Instabilities and the roton spectrum of a quasi-lD Bose-Einstein condensed gas with 

dipole-dipole interactions 



Stefano Giovanazzi^ and Duncan H. J. O'DclP 

'School of Physics & Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, Scotland 
^Dept of Physics & Astronomy, University of Sussex, Falmer, Brighton BNl 9QH, England 

(Dated: February 2, 2008) 

We point out the possibility of having a roton- type excitation spectrum in a quasi- ID Bose- 
Einstein condensate with dipole-dipole interactions. Normally such a system is quite unstable due 
to the attractive portion of the dipolar interaction. However, by reversing the sign of the dipolar 
interaction using either a rotating magnetic field or a laser with circular polarization, a stable cigar- 
shaped configuration can be achieved whose spectrum contains a 'roton' minimum analogous to that 
found in helium 11. Dipolar gases also offer the exciting prospect to tune the depth of this 'roton' 
minimum by directly controlling the interparticle interaction strength. When the minimum touches 
the zero-energy axis the system is once again unstable, possibly to the formation of a density wave. 

PACS numbers: 



I. INTRODUCTION 

Ultra-cold gases of atoms which interact with long- 
range anisotropic dipole-dipole interactions have been 
the subject of a number of theoretical analyses over the 
past few years Illllllllllllllllll|l3. The very re- 
cent advent of Bose-Einstein condensation of molecules 
has renewed interest in the investigation of quantum 
gases with dipolar interactions since molecules can po- 
tentially possess large dipole moments. One of the novel 
properties that has been predicted for gaseous Bose- 
Einstein condensates (BECs) with dipole-dipole interac- 
tions is a 'roton' minimum in the excitation spectrum 
[T3. , a feature which is absent in the usual case of 
repulsive short-range s-wave interactions, and seems to 
originate in the long-range and partially attractive na- 
ture of dipolar interactions. 

The 'roton' minimum can be interpreted with the help 
of Feynman's formula for the dispersion relation for ex- 
citations of energy E and momentum fik of a Bosc liquid 
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where m is the atomic mass. This remarkable formula 
. ^ ' relates the excitation spectrum to the two-particle cor- 
, relation properties as encapsulated in the static struc- 
' ture factor S{k), which is the Fourier transform of the 
5^ ! PE^ir correlation function. A peak in S{k) due to strong 
two-particle correlations can lead to a minimum in the 
energy spectrum. The most famous example of a roton 
minimum (and the system for which the name was first 
coined, by L.D. Landau) is helium II where the roton 
minimum occurs at wavelengths coinciding with the av- 
erage interparticle separation. The roton minimum in 
helium II has often been regarded as indicating that the 
superfluid is very close to crystallization 15, 16, 17] since 
the ordering of atoms on a crystal lattice would lead to 
a corresponding peak in S{k). However, many questions 



concerning this superfluid to density-wave transition re- 
main difficult to answer, essentially because liquid helium 
is a strongly correlated system. Our aim in this paper 
is to explore an analogous 'roton' spectrum in a quasi- 
lD BEG with dipole-dipole interactions. The tunability 
of the various parameters (interaction strength, density, 
effective dimensionality etc) controlling these relatively 
simple quantum gases gives the exciting possibility to 
vary the depth of the roton minimum and hence explore 
the onset of this zero-temperature phase transition. 



II. DIPOLE-DIPOLE INTERACTIONS IN A 
BOSE-EINSTEIN CONDENSATE 

The long-range part of the interaction between two 
dipoles separated by r, and aligned by an external field 
along a unit vector e, is given by 
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Dipoles induced by an electric field E = Ee have a cou- 
pling Cdd — E'^o? I^Q, where a. is the static polarizability, 
and eo the permittivity of free space. For atoms with a 
magnetic dipole moment dm aligned by a magnetic field 
B = Be, one has Cdd = Moc'mi where /io is the per- 
meability of free space. The isotropic short-range part 
the interatomic interactions are modelled by a pseudo- 
potential, C/s(r) = (47ras?i^/m)(5(r) = 5<5(r), where g in- 
corporates the quantum aspects of low-energy scattering 
via the s-wave scattering length, ag. In order to apply 
Bogoliubov (perturbation) theory to a BEG one gener- 
ally needs to use such an effective- or pseudo-potential 
rather than the bare interaction jl^j . In what follows we 
shall assume, however, that the dipole-dipole interaction 
can be treated within the Born approximation, which 
corresponds to the Fourier transform of the bare poten- 
tial |2Jl. The Born approximation has been shown by Yi 
and You to work reasonably well for the long-range 



2 



dipole-dipole interaction, but see Derevianko's work jg 
for corrections, which can be substantial p]. 

The basic quantity that then enters the Bogoliubov 
theory is the Fourier transform of the total effective in- 
teraction, f/tot(k) = / d^rexp[-ik • r]{gS{r) + ?7dd(r)}, 
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where Edd is a dimensionless measure of the strength of 
the dipolar interactions relative to the s-wave scattering 
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The definition includes a factor of 3 because for a homo- 
geneous systeinfno confining trap) one finds instabilities 
when Edd This can be seen directly from 

the Bogoliubov dispersion for phonon-like density pertur- 
bations in a homogeneous dipolar BEC 
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where n is the density. This dispersion relation (jS)) can 
become imaginary when edd > 1, indicating an instabil- 
ity. The interactions enter the dispersion in the second 
term under the square root via the effective potential Q . 
The dispersion Q relation has an angular dependence {0 
is the angle between the momentum of the phonon and 
the external polarizing field) which further illustrates the 
richness of dipolar systems in comparison to the non- 
dipolar case. 

Santos et al considered the case of a quasi-2D pan- 
cafce-shaped BEC, where the dipoles are aligned (by an 
external static field) along the symmetry axis of the trap. 
They predicted that the analogous expression to © gives 
a roton minimum in the excitation spectrum. In this pa- 
per, however, we want to see if a roton minimum can be 
achieved in the quasi- ID case where the condensate is 
cigar shaped, being tightly confined in the radial (x — y) 
plane, say, and very elongated along the axial (z) di- 
rection. At first sight this seems unlikely because the 
anisotropy of the dipole-dipole interaction means that the 
quasi-lD and -2D cases are very different. Two dipoles 
lying side- by-side in the radial plane are repulsive, whilst 
two lying end-to-end along the axial direction are attrac- 
tive so, broadly speaking, a pancake shaped BEC tends 
to be more stable since the mean-field energy is more pos- 
itive (repulsive) than that of a cigar-shaped BEC. Indeed, 
an infinite quasi- ID BEC with only dipole-dipole inter- 
actions (no repulsive s-wave short range interactions) is 
unstable to density perturbations, as we shall see. There- 
fore we shall employ a special trick and reverse the sign 
of the dipole-dipole interaction (see also Q). Two meth- 
ods for achieving the reversed dipole-dipole interaction 
will be reviewed in the next section. With the sign of 
the dipole-dipole interactions reversed, two dipoles lying 



end-to-end along the axial direction are now repulsive 
and in this way an infinite quasi- ID dipolar BEC can be 
stabilized somewhat against collapse induced by density 
perturbations. 

In our previous work on the possibility of a roton spec- 
trum in a BEC ^3 we took the case of laser-induced 
dipole-dipole interactions , which on top of the static 
interactions of Equation |(2Jl, also have a very long-range 
component (l/r^, or even 1/r, depending on the polar- 
ization direction) arising from retardation effects due to 
the finite wavelength of the laser radiation. In princi- 
ple this is a very general scheme since one can bene- 
fit from the large polarizability of atoms/molecules close 
to an electronic transition. The disadvantage of laser- 
induced dipole-dipole interactions, in atomic systems at 
least, is the unavoidable presence of spontaneous emis- 
sion which heats the sample. In molecular systems one 
may be able to circumvent spontaneous emission prob- 
lems somewhat by using transitions in the infra-red or 
even microwave regions where spontaneous emission rates 
are low and the long wavelength of the radiation allows 
one to have many particles in an 'interaction volume' 
(given by Af^gj,^.). In j23 we discussed how a BEC with 
light-induced dipole-dipole interactions is unstable to the 
formation of a density wave. The periodicity of the mod- 
ulation is determined by the wavelength of the laser light 
and in one interpretation can be understood in terms 
of the back-scattered light interfering with the incident 
light to form an optical grating which the atoms then 
minimize their energy in by trying to sit predominantly 
in the wells. In this paper we are interested in the case of 
static dipolc-dipole interactions. Although at first sight 
there are many similarities between the two cases, in the 
light-induced case the optical wavelength provides the 
fundamental (classical) scale in the problem whereas here 
the parameters we must form a length scale from are Cdd, 
n, m and %. 



III. REVERSING THE SIGN OF THE 
DIPOLE-DIPOLE INTERACTION 

i) In the case of magnetic dipoles the sign of the dipole- 
dipole interaction jSJ can be reversed by rapidly rotating 
a component of the external magnetic field in the radial 
plane 5], similarly to a well known technique from solid 
state NMR. The magnetic field should have the form 

B(t) = B {cos (/«z -I- sin 4> [cos {Vlt) x + sin (fit) y] } (6) 

which is a combination of a static magnetic field di- 
rected along z, and a field rotating at frequency 17 in 
the perpendicular radial plane. The frequency should 
be such that the dipoles practically do not move dur- 
ing the time but their moments adiabatically follow 
the external field B(t). This corresponds to the condi- 
tion WLarmor ^ ^ w^, , whcrc Wr, are the radial 
and axial trapping frequencies, respectively. Under these 
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conditions the particles experience a time-averaged in- 
teraction whose form is exactly the same as Equation (j^J 
except for being multiplied by (S cos^ — l) /2. By vary- 
ing the angle (p this factor can be continuously changed 
between -1/2 to 1, giving the possibility to reverse the 
sign of the interaction. At the so-called 'magic angle' 
of (i>M = 54.7° the dipolar interaction averages to zero. 
This technique can be adapted to the case of heteronu- 
clear molecules polarized by a rotating electric field, ii) 
For dipoles induced by the electric field of a very long 
wavelength laser propagating along the z-direction, then 
one obtains a reversed dipolar interaction when the laser 
is circularly polarized. In the static limit the laser in- 
duced dipole-dipole interaction gives 



Cdd - la" Keel) 



(7) 



where / is the laser intensity and c is the speed of light. 
For a laser with wavevector q = iuj/c, then application 

of the identity e*^^^ (q)e^-'^^ (q) = \[{6ij - g^g^) ± ieyfegfc], 
where -f (— ) corresponds to left (right) circular polariza- 
tions, leads to the interaction 



J rcirc 



^'■^ 2 47rce2 



l/a2(0) (l-3cos2 



(8) 



where is angle the interparticle separation vector r 
makes with the z-axis. We see that the sign of the inter- 
action has been reversed and multiplied by one half, cor- 
responding to the maximal change possible in the mag- 
netic case. In order to keep the notation as simple as 
possible, whenever reversed dipolar interactions are used 
in the rest of this paper Cm will be taken to include the 
multiplying pre- factors (such as 1/2). 
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IV. QUASI-ID ANSATZ 



The dipole-dipole energy functional of the gas is given 
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where n(r) and n(k) are the number density and its 
Fourier transform (F.T.), respectively. The F.T. of the 
dipole-dipole potential is given in 10). 

Assuming an harmonic radial trapping Vtrap = 
(l/2)m[j^(a:^ -|-y^) with huox larger or equal to the mean- 
field energy, we make a Gaussian ansatz for the radial 
part of the number density profile 

n{v) = N{i,wl)-^nKz)e^j>[~{x^+y^)/wl] (11) 



H(k) = Nn^{K)cM-{kl + kl)wl/A] 



(12) 



where Wx is the characteristic radial width and N is the 
total number of atoms in the sample. The axial density 
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FIG. 1: Solid curve: effective ID dipolar potential Ul^{v), in 
momentum space in units of Cdd/(27rTO^). Dot-dashed curve: 
a simple approximation (as given in Eq. HIGH ). 



profile n^{z) is normalized to one, and its F.T. is n^{kz). 
Upon inserting the Gaussian ansatz into the dipolar en- 
ergy functional, and taking the polarization direction of 
the dipoles to be in the axial (z) direction the Hm re- 
duces to 

Hm = {N^m j j dzdz' n'^{z)W^K^. - z')nKz%li) 
= {N^/At:) I Ak,W'{k,)nK~k,)UlKkz) (14) 



which defines U^^{kz) as the F.T. of the effective one- 
dimensional dipole-dipole potential. One finds that 
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where = dt exp[— t]/t is the exponential inte- 

gral [23|. In an analogous fashion, the effective ID s- 
wave contact interaction in momentum space is given by 
Ul = 2h'^as/{mwl). 

In order to give a qualitative impression of what the 
effective ID potential looks like in coordinate space we 
(10) note that an analytically simple approximation to U^^ is 



given by (see Fig. ^ 
ULik.) 



dd 



klwl/2 + 1 
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Upon transforming back into coordinate space, the 
l/(fc^w^/2 -I- 1) term gives an exponential contribution 
to the effective coordinate-space ID potential U^^^{z) = 
1/(2^) / dk, exp [ik,z\ U,l,{k,), i.e. 
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V. BOGOLIUBOV EXCITATION SPECTRUM 
FOR A QUASI-ID DIPOLAR BEC 

The Bogoliubov dispersion relation for excitations in a 
homogeneous nearly ideal BEC is 
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where for a quasi-lD system, with a Gaussian radial pro- 
file of peak density n(0), one has = 'nn{{))w'^U^^^/m. 
For the linear parts of the spectrum c plays the role of 
the speed of sound. It is useful to write the Bogoliubov 
relation in terms of the following dimensionless quantities 
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kz is the wavenumber scaled by the inverse radial size, 
Wx, of the condensate radial wave function ansatz. The 
Bogoliubov dispersion relation becomes 

El^fi (^3edd ^kl eMkl]Ei[kl] ^^ + 1^11 + fc^ 

(22) 

When £dd > 1 one finds that the Bogoliubov energy be- 
comes imaginary for a range of momenta, see Figure |21 
indicating an instability due to density perturbations. If 
there are no repulsive s-wave interactions (sdd ^ oo) 
to counterbalance the attractive portion of the dipolar 
interactions then there is an instability even for vanish- 
ingly small dipole-dipole coupling strength in this infinite 
ID case. In an actual experiment the trapping potential 
along the axial direction can stabilize the BEC some- 
what, in an analogous fashion to the case of a negative 
s-wave scattering length p^ . 

However, as pointed out earlier, the stability proper- 
ties of the condensate can be dramatically changed by 
reversing the sign of the dipolar interaction. Figure |3| 
shows the Bogoliubov excitation spectrum in the case of 
negative £dd- In contrast to liquid helium, where one 
only has control over macroscopic thermodynamic vari- 
ables such as temperature and pressure, in an atomic 
or molecular BEC considerable control can be exerted 
over microscopic quantities such as the interparticle in- 
teractions. A well-known example is tuning the s-wave 
scattering length using a Feshbach resonance 23]. Simi- 
larly, in a system of atoms or molecule with a permanent 
magnetic or electric dipole, one can envisage an experi- 
ment at a fixed value of p. and adiabatically changing the 
dispersion relation by adjusting Sdd via the angle of the 
rotating external field. Alternatively, for dipoles induced 
by an electric field one can imagine slowly changing the 
field strength. Within the current model we find that 
as Sdd is increased in magnitude below zero (made more 



FIG. 2: Bogoliubov dispersion relation for axial excitations in 
a quasi- ID BEC with p, — 0.1. Dashed curve: s-wave contact 
interactions only {sdd = 0). Solid curve: dipole-dipole plus 
s-wave contact interactions, £dd = 2. For < 0.19 the Bo- 
goliubov energy for the dipolar case is imaginary. In general, 
whenever £dd > 1 the spectrum is imaginary at long wave- 
lengths, indicating an instability to density perturbations. 
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FIG. 3: 'Roton' minimum in the Bogoliubov dispersion re- 
lation for axial excitations in a quasi- ID dipolar BEC with 
fl = 0.1. Dashed curve: s-wave contact interactions only 
(sdd ~ 0). Solid curve: reversed dipole-dipole plus s-wave 
contact interactions, edd ~ —13. 



negative), first an inflexion point appears in the disper- 
sion relation and then a fully fledged roton minimum as 
shown in Figure |3| As |edd| is increased yet further the 
roton minimum deepens and eventually touches the zero- 
energy axis, at which point, following the earlier discus- 
sion of helium, we might expect an instability towards 
a density-wave. The value of the wavenumber when the 
roton touches the zero-energy axis is found to be 
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VM2(2edd + l)2 + 4/l(edd-l). (23) 




FIG. 4: The critical value of l/£dd(= 3g/Cdd) when the roton 
minimum touches the zero-energy axis plotted as a function 



of /i 



9"(0) 



^ in an infinite quasi-lD dipolar BEG. 



The particular values of the dimensionless chemical po- 
tential ft, which is representative of the density, and 
the dipolar interaction strength, as represented by £ddj 
at which this instability occurs obey the transcendental 
equation 

3/ledd(fcf*)'exp[(fcf = /i(edd-l)-(fcf*)'- 

(24) 

By substituting Equation l|23() into H24|l . one can numer- 
ically obtain the relationship between Edd and p, at the 
instability. The results are shown in Figure 01 

In the case of pure dipolar interactions (no s-wave) the 
point at which the roton touches the zero-energy axis can 
be reduced to a single number 



n(0)C, 
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In this limiting case the wavelength associated with the 
roton minimum is given by Aroton = VSnh/ \J TO7i(0)Cddj 
which sets the scale for the density-wave that we expect 
to form. 



VI. EFFECT OF AXIAL TRAPPING UPON 
THE ROTON INSTABILITY 

Up till now we have discussed the idealized case of 
an infinite quasi-lD cylindrical BEC. We now want to 
move closer to the possible experimental situation where 
the trap is in fact 3D, but is still assumed to be highly 
elongated and have cylindrical symmetry. The ansatz 
for the condensate density may now be taken as a 3D 
Gaussian 



rt(r) = 
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Using this ansatz, the contribution of the dipole-dipole 
interactions to the total energy is found to be given by 



FIG. 5: Dependence of the radial size Wx of a cigar-shaped 
dipolar condensate in a 3D trap upon the system parameters. 



the expression 
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where k — Wxjwz is the aspect ratio of the trapped BEC 
and 



, 1 + 2^2 3K2tanh"^ Vl - 
/W = ^-F ^-21172 • (28) 
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Assuming that the large extension of the trap in the axial 
direction means that the system is in the Thomas-Fermi 
limit for that direction (i.e. the zero-point kinetic en- 
ergy due to the axial trapping can be ignored relative to 
the trap potential and interaction energies — see |1C| for 
a discussion of the Thomas-Fermi limit for dipolar gases) 
then the total energy functional, scaled by the transverse 
trapping energy, is approximately 
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— (Edd - 1) ^ 

Q 



(29) 

where A = lOzI^x is the aspect ratio of the trap, Wx — 
Wx/lx is radial size of the BEC scaled by the oscillator 
length of the trap — ^/fij (mujx), and Wz = Wz/lz is the 
corresponding axial quantity. The first terms in brack- 
ets on the rhs are the radial energy due to the trap and 
zero-point kinetic energy, respectively. The next term is 
the axial trapping energy and the last term is the total 
interaction energy, being the sum of the s-wave and dipo- 
lar contributions. In order to obtain this simple expres- 
sion for -fftot we have expanded /(k) for small n and re- 
tained only the first term. We have also used the quantity 
T] = Ntts/lz which is a key interaction strength parame- 
ter known from regular BECs with contact interactions 
'24^1: when 77 is large the system is in the Thomas- Fermi 
regime. 

Minimising i?tot with respect to variations of the con- 
densate radii gives equations for the equilibrium values 
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FIG. 6: The critical value of l/sdd when the roton minimum 
touches the zero-energy axis plotted as a function of the chem- 
ical potential /i in a very elongated 3D trap. 



of these quantities. One finds 
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Solving these equations gives the condensate radii as a 
function of the various parameters. Figure |31 illustrates 
the dependence of the radius on the parameter combi- 
nation |?7(£dd — 1)1%/^. From Figure El one sees that, all 
other things being equal, the radial size increases as the 
dipolar interaction strength is increased: the cigar be- 
comes fatter. It is useful to note that Equations (|30|l and 
(|31|l allow us to write the total energy functional and the 
actual chemical potential of the BEG in the trap (defined 



as /i = [iJkin 

of the radius 



i/trap + 2ifs + 2iJdd]/iV) solely in terms 

(32) 
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The solutions of equations H3U|) and (|31|l minimise 
the total energy functional in the 2D parameter space 
spanned by the radii w^ and Wz for particular fixed 
values of A, rj and edd- Providing the minimum is a 
global one they correspond to the lowest energy solutions 
one can obtain from scaling variations. Perturbations 
about these equilibrium values should result in stable os- 
cillations which physically correspond to monopole and 
quadrupole shape oscillations. If, for new values of A, 77 
and Edd) the minimum in this {wx,Wz) parameter space 
becomes only a local minimum, or a saddle, then the sys- 
tem is now either only metastable, or unstable, respec- 
tively, to scaling perturbations. The roton instability is in 
principle, however, an independent type of instability due 



essentially to local density perturbations (phonons) in the 
same spirit as the instability we first noted in Equation 
© for a homogeneous system and developed in Section 
Ivl for an infinite cylinder. Local density perturbations 
arc not captured by the simple {wx, Wz) parameter space 
and so must in general be considered in addition to in- 
stabilities arising from scaling perturbations. Scaling in- 
stabilities in trapped dipol ar g ases have been extensively 
discussed before N, a a UJ. Our main purpose here is 
to highlight the local density instability connected with 
the roton minimum touching the zero energy axis. 

We therefore now wish to calculate how the roton in- 
stability point (defined by the roton minimum touching 
the zero-energy axis) depends upon edd and the value of 
the chemical potential in a cigar-shaped 3D trap. In Fig- 
ure 01 we indicated how Edd depends on Ji at the instabil- 
ity. For the case of a 3D trap /2 is not the best quantity to 
consider since it depends upon the density which is itself 
a non-trivial function of Edd via the Equations H3UH31|I . 
In fact one can show 



^^ = 2V2- 

a/tT Wy, £ 



wz. 



dd 



(34) 



For the 3D trap we therefore prefer to illustrate the criti- 
cal value of £dd in terms of the true chemical potential in 
the trap /i, as given by Equation l|33() . The result is shown 
in Figure In the limit that 1/edd ^ we have the case 
of a purely dipolar BEG with no s-wave interactions. We 
find that the limiting value of the chemical potential at 
this point is hKNTilOx) = 1.043 and the corresponding 
value of the condensate radius is Wx = 1.093. These num- 
bers should only be taken as a qualitative guide rather 
than quantitatively accurate since the BEG is no longer 
one dimensional when ^/ {Nfiojx) ~ 1, and the effective 
ID potential approximation will begin to break down. 
In particular, one should include radial excitations in or- 
der to give a consistent treatment. This is beyond the 
scope of this largely illustrative paper, but for a pancake- 
shaped dipolar BEG such a treatment can be found in 
|13|. Nevertheless, the very simple treatment given here 
should give the general picture of what to expect. 



VII. CONCLUSION AND OUTLOOK 

A infinite quasi- ID dipolar BEG gas would normally 
be unstable to local density fluctuations if e^d > 1 on 
account of the attractive nature of the dipolar interac- 
tions along the axis of the trap. However, by reversing 
the sign of the dipolar coupling using either a rotating 
polarizing field for permanent dipoles or a very long 
wavelength circularly polarized laser beam for electrically 
induced dipoles, one can obtain a stable quasi-lD system. 
This configuration has the possibility of a 'roton' mini- 
mum its long-wavelength axial excitation spectrum. This 
minimum is tunable via parameters such as the dipolar 
interaction strength and the density. As is suspected to 
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be the case in helium II, we speculate that the roton min- 
imum could be the precursor of a transition to a density 
wave. The tunability of dipolar gases means that this 
phase transition could be explored experimentally when 
quantum gases with significant dipole-dipole interactions 
are realized. One of the important features of the ID sys- 
tem discussed in this paper is that the change in the order 
parameter between a superfluid and a density wave does 
not seem to involve any fundamental change of symme- 



try and so the phase transition could be a smooth one of 
2nd order. For a pancake or fully 3D system the change 
in symmetry is more dramatic and these transitions may 
not be smooth. We believe considerations such as these 
make dipolar quantum gases very worthwhile systems to 
study. We acknowledge financial support from the EP- 
SRC (UK) and the Marie Curie Programme of the Euro- 
pean Commission. 
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